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It is shown that varying speed of light cosmology follows from a string-inspired minimal length 
uncertainty relation. Due to the reduction of the available phase space volume per quantum mode 
at short wavelengths, the equation of state of ultrarelativistic particles stiffens at very high densities. 
This causes a stronger than usual deceleration of the scale factor which competes with a higher than 
usual propagation speed of the particles. Various measures for the effective propagation speed are 
analyzed: the group and phase velocity in the high energy tail, the thermal average of the group 
and phase velocity, and the speed of sound. Of these three groups, only the first provides a possible 
solution to the cosmological horizon problem. 
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I. INTRODUCTION 

Varying speed of light (VSL) cosmology was conceived 
as an alternative solution to the cosmological horizon 
problem ^, IJ]. In contrast with inflation, where 

the comoving horizon scale shrinks as a consequence of 
the accelerating growth of the scale factor, VSL keeps the 
scale factor dynamics of the standard hot big bang un- 
changed but assumes that the speed of light, c{t), is a de- 
creasing function of of cosmic time t. The horizon prob- 
lem can be solved provided that c{t) drops sufficiently 
steeply for some time. In the original formulation, c{t) 
was treated as a free parameter of the theory. 

In Ref . 1^ , it was noted that VSL can effectively be ob- 
tained from the thermodynamical properties of theories 
with nonlinear dispersion relations which, in turn, were 
motivated in |^ by noncommutative geometry. While 
there are indications from string/M-theory that noncom- 
mutativity of spatial coordinate operators may be rele- 
vant near the string scale, we are lacking a simple, con- 
crete realization of noncommutative geometry that can 
be used in the cosmological context. Hence one is forced 
to reduce its effects to a modified dispersion relation. 

A very similar yet more specific approach is taken in 
this work. The starting point is an observation that 
was made in numerous studies of quantum gravity and 
string theory: under rather general assumptions, one can 
show that there exists a minimal length scale that can be 
probed by experiments whereas the resolution of shorter 
distances is prohibited by quantum gravitational effects 
(e.g., [|[ 0, H, Examined from the point of view of 
the low energy effective theory, it appears that quantum 
gravity modifies the Heisenberg uncertainty relation to 

AxAp > ^ (1 + + ■■■) (1) 

which gives rise to a minimal short distance uncertainty 
Axinin = even for infinite Ap. For practical pur- 

poses, Axniin can be identified with either the string scale 
or the Planck length, where the former could be a few or- 
ders of magnitudes larger than the latter. The correction 



may of course originate from complicated dynamics of 
the underlying fundamental theory. On the other hand, 
it can be modeled by equivalent corrections to the canon- 
ical commutation relation: 



[x,p] = z;i(i-f /3p2 + ...) 



(2) 



as discussed in Indeed, the uncertainty relation 

above belongs to class of only very few types of short- 
distance structures of space-time that are admitted under 
very general assumptions JlTf . Hilbert space representa- 
tions of Eq. H were employed, for instance, for regulariz- 
ing field theory and, more recently, for analyzing the 
impact of short distance uncertainty on the predictions 
of inflation pi |1. 



The flnite spatial localization of particles also changes 
their thermodynamic behavior at very high density. This 
has already been observed in |l^ where some cosmo- 
logical consequences have been pointed out. However, a 
complete analysis of this theory in the context of VSL 
cosmology has not been done and shall be provided in 
this work. The most relevant features of the minimal 
length uncertainty theory are summarized in Sec. ||, and 
are used to derive the expressions for the energy density 
and pressure of a radiation fluid in Sec. III. In Sec. pv| , 
the cosmological implications are discussed and a further 
outlook is presented in Sec. 0. 



II. QUANTUM THEORY WITH MINIMAL 
SHORT DISTANCE UNCERTAINTY 



An explicit example for the choice of commutation re- 
lations in three spatial dimensions that maintains trans- 
lation and rotation invariance is given in UM- 



[Xi,pj] = ih 



(l-f 2/3p2)i/2 _ 1 
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Demanding translation and rotation invariance, this 
choice is unique to first order in /3 ^. It does, however, vi- 
olate Lorentz invariance and thereby specify a preferred 
frame which, in the present context, coincides with the 
cosmological rest frame. 

Owing to the finite short distance uncertainty, the 
usual Hilbert space representation in terms of positions 
eigenstates is no longer available (in contrast with the 
momentum representation \p), which still is). An alter- 
native representation is found by introducing the trans- 
lators with eigenstates \k) obeying T^Ik) — Ki/ih\K), 



i = 1 . . . 3, and < 2//3 
mcntum space is given by 

{k\p) = (5 (pi - 



The projection onto mo- 



(6) 



In terms of ^{n) — {k\^P), the operator representations 
and the scalar product become: 



1 - /3k2/2 



(V'l|^2 



d^K ■0i('«)'02(k) 



(7) 

(8) 

(9) 



^<2//3 



The same representation was employed in . 

As demonstrated in , it is now possible to construct 
the closest analogue to the position representation by re- 
placing the usual position eigenstates with states of max- 
imal localization around a given position, |a;™'). Their k- 
space representation can be constructed by applying the 
translation operator e^"^ to the maximally localized field 
around the origin, and therefore it varies as ~ Q-txK./h_ 
Using (j^), one finds the quasi-position representation of 
the plane wave with momentum p: 



{p\x 



ml\ 



N{p^) exp 



-ixpyjl + 2^p2 _ I 



(10) 



where N(p'^) is a normalization coefficient. It is now clear 
that k/H acts as a substitute of the usual wavenumber k, 
with the crucial difference that it only allows a minimum 
wavelength Amin = 7r?i(2/3)^/2 for p — > oo, as required at 
the onset. This will become important for the computa- 
tion of the thermodynamic variables in the next section. 



III. EQUATION OF STATE OF A RADIATION 
FLUID 

In order to derive the energy density and pressure of an 
ideal gas of photons (with straightforward generalization 



^ Note that in |l5|, a different choice for the commutation relations 
was made which is noncommutative in the spatial coordinates 
and breaks translation invariance. 



to other ultrarelativistic particles), the grand canonical 
formalism can be used (e.g. ||l^). The grand potential 
is defined as 



n = -k^T InZ , 
where Z is the grand partition function. 



(11) 



(12) 



H is the Hamiltonian, N is the number operator, fee is 
Boltzmann's constant, and T is the temperature. Evalu- 
ating the trace in the number representation for bosons, 
one finds as usual 



1 



-(E''-ii)/kuT 



(13) 



The summation is over all allowed momentum states 
and = c|p'| is the energy of a relativistic particle in 
this state. 

At this point, one usually transforms the sum into an 
integral over the three-dimensional wavenumber k by par- 
titioning space into boxes of volume V — and demand- 
ing L-periodicity of plane waves with momentum p = hk 
in the position representation. In the present framework, 
a position representation is unavailable. Instead, the con- 
tinuum limit is now most conveniently taken in K-space: 
requiring the periodicity of e'^'*/'' yields n ~ 2TThl/L, 
where I € M labels the discrete quantum state and 
Af = {/ e Z"|P < (i/Amin)^}. For sufiiciently larse 
V, the energy levels are nearly continuous and Eq. ( |l3| ) 
can be written as 



gVk^T 



d'^K In 1 



-[E{K)-f^]/kBT 



(14) 



K^<2/I3 



with E{k) = c|p(k)| and g = 2 for photons. After 
transforming the integral to energy space and writing 
it in terms of the dimensionless variables e — E/Ep, 
jji — ^/Ep, and r — kBT/Ep, where Ep = c(3~^/'^ is the 
characteristic energy scale where quantum gravitational 
effects become relevant, one obtains 



Vt = 



AirgVE^T 
{2TThcf 



OO 

y"deJ(e)e2 ln(l-e-['-'^l/") . (15) 



Here, J is the Jacobian determinant of the transforma- 
tion Ki ^ Pi = Ki{l — i3k^/2)~^, Pi e M, given in terms 
of e by 



J(e) 



2-l-e^ 



2 + 3e2 

VTT2? 



(16) 



J(e) essentially contains all of the information about the 
modification of the phase space volume that results from 
the short distance cutoff. 

One can now derive the thermodynamic variables from 
the grand potential in the standard way. Using pir) = 



3 



9t--i T~^rt for the energy density, P{t) = —dv ^ for 
the pressure, N{t) = -~V~^ for the number density 
and setting /2 = for ultrarelativistic particles, one finds: 



P{r) 
P{t) 
N{t) 



(27rfic)3 



de J(e) 



(17) 



[ de J(e) / ^ 



(19) 



It is useful to study the asymptotic behavior of the 
integrals for very small and very large r (e.g., |18 ). In 
the case of p(t), rescaling the integrand by substituting 
u = e/r shows that the dominant contribution comes 
from u ~ 0{1) for r ^ 0. One can therefore write J = 
l+0{u^T^) and find the usual result p{t) ^ for r — > 0, 
as expected. For large r, we first note that J(e) drops 
as e~'* for large e. Hence, for e = 0(1) the integrand 
is 0(t) and so is the integral, whereas for e = 0(t), 
the integrand is 0{t~^) and the integral is 0(1). The 
dominant contribution therefore comes from e = 0(1) so 
that the exponential can be approximated by 1 + e/r + 
0(e^r~^), yielding 



J de J(e 



(20) 



where the integral evaluates to 2^/^/3. 

In the case of P{t), partial integration yields 



P{t) 



3(27rftc)3 



deG(e). 



(21) 



where G(e) is defined as 



G(e) = deJ(e), 



= Vl + 2e2 J(e) 



(22) 



For small r, G(e) = 1 — 0{e^) behaves like J(e), and one 
obtains the expected result P — p/3. On the other hand, 
G ^ e~3 drops more slowly for t ^ oo than J(e). It 
is only for this reason that the ratio of pressure and en- 
ergy density, and hence the dynamics of the scale factor, 
differs from its usual behavior at high temperatures, as 
will be demonstrated in Sec. 
the integrand is now 0(1) 



IV 



and t 



Indeed, for e = 0(t), 
le integral is 0(t), con- 
In 



tributing equally strongly to the integral as e = 0(1) 
this good approximation is given by 



0(r 

Pr^oo(T) - J 



de- 



0(1) 

t{A + B log r) 



(23) 



The cosmological evolution of p is determined, via the 
energy conservation equation, by the equation of state 
parameter w = P/p. For small r, one recovers the 
usual equation of state, w = 1/3, whereas r — > oo gives 
Wr^ooip) — C + D\og Pt-^oo- The same form of w was 
discovered and analyzed in hence one can expect a 
similar result for the cosmological evolution. As shown in 
Sec. IV, this is indeed true. The result for w{p) was con- 
firmed numerically by inverting p{t) in order to evaluate 
P{t). a good high-density fit is provided by G ~ —1.2 
and D ~ 2.3. 

Finally, using the same arguments as for p{t), the ex- 
pression for N{t) yields the standard result for r — 
while T OO gives 



A^r^oc.(T) 



{2'KhcY 



OO 

j deJ(e)e 



27r g^l 
(27rftc)3 



(24) 



All of the preliminaries are now in place to examine the 
cosmological impact of the modified equation of state, 
which will proceed in close analogy with Ref . [H . 



IV. COSMOLOGICAL IMPLICATIONS 

The evolution of p{t) and the scale factor a{t) in the 
very early universe are determined by the Friedmann and 
energy conservation equations with negligible curvature 
and cosmological constant terms, given by 



3Af2j 



P 



--{l + w{p)) , 
a 



(25) 
(26) 



where a dot corresponds to a derivative with respect to 
cosmic time and Mpi = (SttG)"^/^ is the reduced Planck 
mass. These equations can only be solved numerically 
(starting with the known low-temperature solution and 
integrating backward) , using the numerical evaluation of 
w{p). However, if w(p) varies sufficiently slowly with p, 
one can approximate the solution for the scale factor by 
the usual one for constant w, with w replaced by w{p) 
and an offset for the origin of the time coordinate to: 



ait) ~ (t - to)2/3(l+-(p)) 



(27) 



Fig. |l| demonstrates that this is a good approximation at 
all times except for the brief period where w transitions 
from 1/3 to the logarithmic high-temperature behavior. 

Eq. ( |27| ) captures the relevant scale factor dynamics 
in our model: at very high density, the deceleration ~ 
d/aH^ is very large, giving rise to a smaller than usual 
Hubble rate H — d/a, whereas the standard scale factor 
evolution is recovered at densities far below the string 
scale. 

The high-density evolution of p as a function of the 
scale factor can be obtained directly by solving Eq. (|2l) 
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FIG. 1: Comparison of the effective time exponent of the 
numerical solution (solid line) and the approximation 2/3(1 + 
?ii(p)) (dashed line) as a function of time, adjusted for a shift 
of the time of the cosmic singularity. 



with w(p) ~ ifr^oo, yielding 
logp^^oo(a) ^Ea 



-ZD 



C + 1 
D 



(28) 



This was also shown in ||] . 

The question whether the cosmological horizon prob- 
lem can be solved in this model hinges crucially on the 
choice for the effective propagation speed of information, 
Ceff. The VSL mechanism works if the comoving horizon 
scale, 



Ceff 

aH 



Ceff 

d 



(29) 



declines for some period of time before turning around 
to its current growing behavior. 

A case can be made that the most suitable measure 
for Coff is given by the speed of sound, Cg, of the radia- 
tion fluid rather than by the group or phase velocity of 
individual photons. It is given by 



dp_ 

dp 



S— const 



dw 
' dp 



(30) 



In the high-density limit, Cg scales like (logp)^/^ ^ a~^/^, 
whereas a ^ ap^^^ ~ aexp(l/2a'^) which clearly wins as 
a — > 0. Consequently, is always a growing function of 
time in this case and the horizon problem remains. 
The group and phase velocities are given by 



Cg(e) 



Cph(e) = 



dE 
E 



(l + 2e2 + VTT2^) (31) 
l + Vl + ^eA . (32) 



As pointed out in , the evaluation of e in these expres- 
sions is somewhat ambiguous. The bulk of all photons 
will populate a region around the peak in the modified 



Planck distribution, which saturates at r ~ 1. After 
saturation, the group and phase velocities of these pho- 
tons will cease growing as a function of increasing den- 
sity. However, the authors of |^ argue that the "fast 
tail" of the Planck distribution may provide the required 
causal contact, and they propose to evaluate Cg and Cph 
at e ~ T ~ ►oo- In this case, Cg ~ p^ and Cph ~ p at 
high densities, so that rh indeed declines as a function 
of time (consistent with |Q). Consequently, this choice 
of effective velocity solves the horizon problem, provided 
that the model assumptions made in Sec. || remain valid 
well into the regime r 3> 1. 

As a final possibility, one can consider a thermal av- 
erage of Cg (or Cph), defined for instance as the particle 
number weighed value: 

(cg,ph)(r) = N-\r)j,^ /d,-M±U_(33) 



{2nhc) 



Kg^ph{e) = Cg,phJ(e) 



(34) 



The asymptotic properties of N{t) were discussed in 
Sec. III. Noting that A'g.ph = 1 — O(e^) for small e one 
finds (Cg_ph)(T) ~ c as required. On the other hand. 
Kg ~ for large e, whereas Kph ^ e^^. The situation 
for the integral over Kg at r ^ oo is therefore analogous 
to that of Pr^ao, yielding the same logarithmic growth. 
The normalization by N{t) gives rise to the same behav- 
ior as found earlier for W'j- — >oQ, i.e. (cg)^- — >oo ^ P- In 
contrast, the integral over ifph at r — s- oo is analogous to 
that of Pt^oo, hence the normalized result for (cph)T^oo 
asymptotes to a constant independent of t. Explicitly, 
one finds (cph)T-^oo — 2c. 

To summarize, the thermally averaged group or phase 
velocities grow at best logarithmically. This is insuffi- 
cient to solve the horizon problem by means of the VSL 
mechanism, as argued above for Cg. 

Fig. ^ shows the size of the comoving horizon as a 
function of time, computed numerically using the "fast 
tail" group velocity and the thermally averaged one. 



V. SUMMARY AND OUTLOOK 

The analysis presented in this work is based on a 
rather simple premise: The spatial separation that can 
be probed by high-energy particles has a finite minimum 
value. Without referring to any specific quantum grav- 
itational effect in particular, it is a very general model 
for the physics of relativistic particles at densities that 
occurred in the very early universe. It is very interesting 
that such a simple modification may provide an alterna- 
tive solution to the horizon problem without resorting to 
the inflationary paradigm. Of course, one must now think 
of possibilities to resolve other quandaries that inflation 
has proved to be so helpful with, such as the flatness and 
relics problems and, perhaps most importantly, the gen- 
eration of scale invariant fluctuations. Some ideas in this 
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FIG. 2: Comoving horizon distance as a function of scale 
factor for the "fast tail" group velocity Cg (solid line) and the 
thermally averaged one, (cg) (dashed line). While the former 
can solve the horizon problem, the latter cannot. At late 
times, both converge onto the radiation dominated solution 



direction have already been put forward in the context 
of VSL § |, and they can be readily generalized to the 
model in this work. 

As always, many open questions remain. One of them, 
which may be answered by a more detailed analysis of 
transport processes in the early universe, is the correct 
choice of the effective communication speed. On a more 
fundamental level, one may question the validity of the 
minimal distance uncertainty principle as early as may be 
required to solve the horizon problem. In other words, 
it is possible (and even likely) that additional, quantum 
gravitational degrees of freedom will be excited in the rel- 
evant regime of p that will drastically change the equation 
of state. Unfortunately, the answer to this question might 
remain elusive until the fundamental theory is known in 
full detail. 
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